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Abstract 

Three-scale homogenization procedure is proposed in this paper to provide estimates 
of the effective thermal conductivities of porous carbon-carbon textile composites. 
On each scale - the level of fiber tow (micro-scale), the level of yarns (meso-scale) 
and the level of laminate (macro-scale) - a two step homogenization procedure based 
on the Mori-Tanaka averaging scheme is adopted. This involves evaluation of the 
effective properties first in the absence of pores. In the next step, an ellipsoidal 
pore is introduced into a new, generally orthotropic, matrix to make provision for 
the presence of crimp voids and transverse and delamination cracks resulting from 
the thermal transformation of a polymeric precursor into the carbon matrix. Other 
sources of imperfections also attributed to the manufacturing processes, including 
non-uniform texture of the reinforcements, are taken into consideration through the 
histograms of inclination angles measured along the fiber tow path together with a 
particular shape of the equivalent ellipsoidal inclusion proposed already in (Skocek 
et al., 2008). The analysis shows that a reasonable agreement of the numerical 
predictions with experimental measurements can be achieved. 

Key words: Carbon-carbon composites, multi-scale analysis, Mori-Tanaka 
method, optimization, porous materials 



1 Introduction 



Since its introduction the Mori-Tanaka (MT) method (Mori and Tanaka, 1973) 



has enjoyed a considerable interest in a variety of engineering applications. These 
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(a) 
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Fig. 1. Color images of a real composite system: (a) Scheme of multi-scale structural model 
(from top - transverse and longitudinal view of fiber tow composite, composite unit cell, 
composite lamina, composite plate), (b) Carbon tow microstructure showing major pores 
and transverse cracks. 



include classical fiber matrix composites (Benveniste 1987; Dvorak and Sejnoha 



1996), natural fiber systems (Hellmich and Ulm, 2005), or even, although to a lesser 



extent, typical ci vil engineering materials suc h as asphalts (Lackner et al, 2005) 



or cement pastes (Smilauer and Bittnar 2006). While generally applied to modest 



geometries, it has also been demonstrated that the MT method may reliably assist 
the engineer to meet the ever growing challenge associated with the analysis of new 



highly complicated material systems such as textile reinforced composites ( Gommers 



et al. 1998 Huysmans et al, 1998). 



Recent studies addressing the behavior of imperfect carbon-carbon (C/C) textile 



composites ( Skocek et al. , 2008 ) further promote the use of the MT method as it 
allows for a direct introduction of various types of imperfections, e.g. in the fiber 
tow path represented for example by histograms of distribution of the fiber-tow ori- 



entation angles (Koskova and Vopicka, 2001). Although carbon-carbon plain weave 



composites belong to a progressive material systems with many applications, a rele- 
vant micromechanics model taking into account most of the geometrical details is still 
missing. Several appealing routes have been already offered to satisfactorily reflect 
commonly observed imperfections both in the woven path and through the laminate 



thickness developed during the manufacturing process (Zeman and Sejnoha 



2004 



Skocek et al. , 2008). The reported works failed, however, to include one of the most 



important features of C/C composites illustrated in Fig. [T]- the porosity, which in 
real systems may exceed 30% of the overall volume. 



A general awareness of the need for incorporating the porous phase in the predic- 
tions of overall response of C/C composites has been manifested in several recent 
works. While all microstructural details were properly identified, the actual analysis 
was limited to either unidirectional fiber composites represented here by individual 



yarns ( 


Tsukrov et al.. 


2005 


Fiat, Tsukrov, Mladenov, Verijenko, Guellali, Schnack 


and Hoffman 


20071 


Fiat, Tsukrov, Mladenov, Guellali, Ermel, Beck and Hoffman, 



2 



2007) or finite element simulations of entire laminate performed in two-dimensional 



(2D) environment only (Tomkova et al., 2008). An extension of this topic taking 
into consideration the characteristic three-dimensional (3D) structure of C/C tex- 
tile composites is presented in this paper. The formulation given here is in the spirit 



of multi-scale analysis discussed in (Tomkova et a/., 2008) combined with the ap- 



plication of the MT method to the prediction of effective elastic properties of C/C 
composite presented in (Skocek et a/., 2008). However, due to the constant impor- 



tance of high-temperature systems with a particular role of C/C composites we 
consider in this paper the subject of effective thermal conductivities. 

The paper is organized as follows. General description of the Mori-Tanaka method 
in the framework of steady state heat conduction problem is outlined in Section |2} 
For an extensive list of references in this subject the reader is referred to (Hatta and 



Taya 1986 Benveniste et al, 1990 Jeong et al, 1998). Two particular issues are 



addressed: solution of the problem of a solitary ellipsoidal inclusion embedded into 



an orthotropic matrix, see e.g. (Chen and Yang 1995), and evaluation of orientation- 



dependent average fields. The ordering of the remaining sections follows the concept 
of the assumed uncoupled multi-scale homogenization approach in which the results 
derived from the homogenization step on a lower scale are used as an input to the 



same analysis performed on the upper scale. Following (Tomkova et al] 2008) three 
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particular scales are examined. The level of fiber tow evident from Fig. Ilj(b) is treated 
in Section |3] while Section |4] examines various geometrical scenarios encountered 
at the level of textile ply, see Fig. [l|^a). Section |5] then provides the estimates of 
the effective thermal conductivities of the laminate and compares those with the 
available experimental measurements. Standard matrix notation is used throughout 
the paper. 



2 Mori-Tanaka method 



As illustrated in Fig. [Tjboth micro- (fiber tow) and meso- (textile ply) scales call for 
treating at minimum a three phase composite medium comprising a homogeneous 
matrix, a certain type of reinforcement and voids. Combining all three phases in 
a single homogenization step appears, however, rather inadequate owing to a con- 
siderable difference in both size and shape between the reinforcements and voids. 
Therefore, a two step homogenization approach is adopted in this paper. In partic- 
ular, the effective properties of a composite aggregate are first found in the absence 
of pores followed by the second homogenization step in which the porous phase is 
introduced into a new homogenized matrix. 

It is therefore sufficient to consider a two-phase composite medium with the hetero- 
geneities (reinforcements or voids) having in general a certain orientation distribu- 
tion. Let this composite be subjected on its outside boundary F to a homogeneous 



temperature boundary condition defined as (Benveniste et al, 1990 Tomkova et al. 
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Fig. 2. Local coordinate system and definition of tlie Euler angles. 



2008) 



(X) = H^X 



on 



(1) 



where H represents the macroscopically uniform temperature gradient vector and 
X are the components of the assigned (fixed, global) Cartesian coordinate system, 
Fig. [2| The volume average of the local constitutive equation written in the local 
a;-coordinate system is then provided by 



(2) 



where q^. is the volume average of the local heat flux in the phase r (r = 1,2 
with r = 1 reserved for the matrix phase) and Xr is the corresponding thermal 
conductivity matrix [Wm~^K~^]. Also note that x = X for the matrix phase. 

Next, in the context of the Mori-Tanaka method, consider a certain auxiliary trans- 
formation problem where a single heterogeneity in an infinite matrix is replaced 
by an equivalent inclusion of the same shape and orientation but having the ma- 
terial properties of the matrix phase. In the Mori-Tanaka mean field theory the 
mutual interaction of individual heterogeneities is taken into account by loading 
such a configuration at infinity by the average temperature gradient in the matrix 
hi, see e.g. (Benveniste, 1987; Benveniste et a/., 1990) for further details. The local 



temperature gradient in the inclusion is then expressed as 



"■2 



(3) 



where h* represents a certain transformation temperature gradient introduced in 
the homogeneous, generally anisotropic, matrix to give the same local fluxes as in 
the composite. Note that the second order tensor S and the vector h* are analogous 
to the Eshelby tensor and transformation strain known from the solution of the 
Eshelby equivalent inclusion problem in elasticity (Eshelby, 1957[ ) (hereafter, we 
shall quote the "Eshelby problem" whenever referring to a general transformation 
inclusion problem). Recall that in the MT scheme the S tensor is a function of 
the matrix material parameters and the shape of the inclusion. Combining Eqs. ^ 
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and (|3j) for the heterogeneity and equivalent inclusion then gives 

92 = -X2^2 = -Xl {hi - {h*r) , 

so that 



{h*r = {xl)-'{xl-X'2)hl 
Next, substitute Eq. ^ into ([s]) to get 

hl = hl + S-{xl)-'{xl-xl)hl 



and finally 



hi 



''•1 ~ ' 2'''l5 



(4) 
(5) 

(6) 
(7) 



where I is the identity matrix and the matrix T2 is referred to as the partial con- 
centration factor. 

Suppose that the heterogeneity possess a certain orientation described by the ori- 
entation distribution function g{(f),6,(3) with (j),6 and (3 being the Euler angles. A 
particular form of g for plain weaved composites is given later in Section |4| In gen- 
eral, following (Jeong et al., 1998), the overall average temperature gradient for a 



two-phase composite with an orientation-dependent inclusion given in the global 
X-coordinate system then attains the form 



{h) = cihi + C2 ((^2) 



(8) 



where is the volume fraction of the phase r, ( ) stands for the volumetric averaging 
and the double brackets (( )) denote averaging over all possible orientations. The 
vector /12 in Eq. ([s]) follows from standard transformation of coordinates so that 

(9) 



X = Qx and ^2 = Q^o = QTnQ^^i = T2^i 



A specific form of the transformation matrix Q consistent with Fig. i M reads 



Q 



cos 4> COS 9 cos 13 — sin (j) sin /? sin 4> cos 9 cos /? + cos (j) sin /? — sin 6 cos /3 
- cos (j) cos 6 sin [3 — sin (j) cos (3 — sin (j) cos 9 sin j3 + cos cos (3 sin 9 sin (3 
cos (b sin 9 sin 6 sin 9 cos 9 



Next, suppose that the local temperature gradient /i2 is expressed in terms of the 
prescribed macroscopically uniform temperature gradient H = {h) as 

h2 = A2/f, (10) 

^ Note that so-called "x2 convention" is used; i.e. a conversion into a new coordinates 
system follows three consecutive steps. First, the rotation of angle (p around the original 
X3 axis is done. Then, the rotation of angle 9 around the new X2 axis is followed by the 
rotation of angle P around the new X3 axis to finish the conversion. 
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(a) 



(b) 



Fig. 3. Homogenization on micro-scale: (a) fiber tow composite, (b) fiber-matrix composite 



where the matrix A2 is termed the concentration factor (Benveniste et al. 1990 



Jeong et al] 1998|) for the heterogeneity. Clearly, the orientation average of /12 then 

(11) 



follows from 

((^2)) = ((T2)) h, 
which, together with Eqs. ([7])-([9]) gives 

{{fi^2)) = {{T,))[cJ+C2{{T2))r'. 



Combining Eqs. (|8j) and (11) further gives 

C,h, = [\-C2{{A2))]H. 



(12) 



(13) 



In analogy to Eq. ([8j), the volume average of the overall heat flux is provided by 

{q) = c,q, + C2 ((^2)) • (14) 
Next, define an overall conductivity matrix x in the fixed coordinate system X and 

XH = CiXlhl + C2 ((X2^2)) • 



write Eq. (14) as 



Introducing Eqs. (10) and (13) into Eq. (15) finally provides 



X = Xi + C2 [((X2A2)) - Xi ((A2))] . 



(15) 
(16) 



Particular forms of the homogenized matrix x be now given for individual 
micromechanics problems. 



3 Micro-scale 



In the first step of the proposed multi-scale homogenization scheme we consider a 
single filament (fiber tow) of a plain weave carbon fabric Hexcel 1/1 bonded to a 
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carbon matrix, Fig. [sj^a). Each filament contains about 6000 carbon fibers T800H 
and significant amount of transverse cracks and voids resulting in a porosity of more 



than 10% (Tomkova and Koskova, 2004 Tomkova et al., 2008). While the matrix 
phase, which essentially corresponds to a glassy carbon, is assumed isotropic the 
carbon fiber possess a transverse isotropy with the value of longitudinal thermal 
conductivity considerably exceeding the one in the transverse direction. The phase 
thermal conductivities are listed in Table [H Considerable difference in the size of 
the two types heterogeneities (fibers and pores) readily suggests a two step homog- 
enization procedure to predict the effective properties of the fiber tow as discussed 
next. 



Table 1 



(TORAYCA 


2008 


Ohlhorst 


1997 



Material 


Thermal conductivity 
[Wm-iR-i] 


Carbon fibers 
Carbon matrix 
Voids filled with air 


(0.35, 0.35, 35) 
6.3 
0.02 



3.1 Effective conductivities of fiber matrix composites 



Fig. |3](b) shows a representative section of the fiber ma trix composite taken from 
the fiber tow in Fig. 3j(a). Based on our previo us studies (Zeman and Sejnoha 



Sejnoha and Zeman 



2002 



Sejnoha et al. 



2004 



2001 



to cite a few) such a composite can be 
quantified as ergodic, statistically homogeneous with a random distribution of fibers 
having the volume fraction of approximately 50%. In the MT scheme the effective 
properties follow from the solution of an auxiliary problem where an infinite cylinder 
of a circular cross-section with semi-axes — > oo, ^2 = ^3 (the xi axis assumed in 
the fiber direction) is embedded into an infinite isotropic matrix. In this particular 



case the effective thermal conductivity matrix given by Eq. (16) simplifies as 



X = Xl + C2 (x2-Xi)A2, 



(17) 



where 



A2 = T2 [cil + C2T2] and T2 = l-Sxi^(xi-X2 



(18) 



An explicit form of S for this particular case of aligned circular fibers in an isotropic 



matrix is available in (Hatta and Taya 1986). 
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Having derived the effective properties of the fiber matrix composite we proceed 
with the second homogenization step to account for the porous phase. As evident 
in Fig. [Sj^a), several distinct shapes of voids can be identified. It is certainly out 
of the question to treat each void separately. Therefore, in the present study, they 
are all combined into a single equivalent inclusion resembling an elliptic cylinder. 
Here, the cylinder is embedded into a transversely isotropic matrix. However, since 
the 5*11 component of S is equal to zero, the solution of an elliptic cylinder in an 
isotropic matrix summarized in (Hatta and Taya 1986) is again applicable. This 



results in the same form of estimate of the effective conductivity matrix x as given 



byEq. (17). 



Nevertheless, there is still one open problem associated with the shape of the el- 
liptical cross-section. Clearly, since the equivalent inclusion represents all possible 
shapes of voids, it can hardly be determined directly from the images of real com- 
posites such as the one in Fig. [sj^a). Instead, to solve this particular problem, we 
exploited the results available from the finite element (FE) simulations carried out 
in (Tomkova et al., 2008). In particular, the optimal aspect ratio ^2/^3 (^1 C)o) of 
the elliptical cross-section was found by matching the effective material properties 
derived from both the MT method and FE solutions. Fig. |4] shows a variation of the 
objective function E 



E 

.i=l 



.1/2 



FEM 

A.ii 



MT 
Am 



(6,^3 



The resulting effective thermal conductivities provided by the proposed two-step ho- 
mogenization scheme are stored in Table |2] (Note that the fiber volume fraction c/jber 
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was estimated from Fig. |3]^b), while the volume fraction of voids Cyoid = Vvoid/Vtow 
stands for the total volume of voids in the fiber tow in Fig. [3]^a)). It also clearly shows 
a significance in properly choosing the shape of the cross-section of the equivalent 
elliptic cylinder. Thereby, to make the analysis more robust, an empirical relation 
between the observed porosity and representative equivalent inclusion is needed. 
This particular topic enjoys our current research interest. 

Table 2 

Effective thermal conductivities of the fiber tow 



Material 


Equivalent inclusion 


Thermal conductivity 
[Wm-iR-i] 


fiber {c fiber) 


void {cyoid) 


Fiber-matrix composite 


oo, 1, 1 (0.55) 




22.09, 2.14, 2.14 


Porous fiber-tow 
(optimal aspect ratio) 




oo, 1, 1 (0.12) 
oo, 1, 10 (0.12) 
oo, 1.6, 10(0.12) 


19.44, 1.54, 1.54 
19.44, 1.02, 1.63 
19.44, 1.12, 1.85 



4 Meso-scale 



While unidirectional fiber matrix composites reviewed in the previous section have 
been of a general interest since some fifty years ago, composite systems with a 
formidable textile texture have received more attention from both academic and 
industrial communities only recently. 

This section examines, at least from the geometrical point of view, the most simple 
representative - a plain weave textile composite. One of the earliest known compu- 



tational models focusing on actual geometry of the textile ply is developed in (Kuhn 



and Charalambides 1999). A three-dimensional view with a typical cross-section are 
plotted in Fig.|5](a)(b). The idealized geometry of this model assumes the centerlines 
of the warp and fill systems of tows to be described by a simple trigonometric form 

c(x) = -sin(^— j. (19) 



Although tempting, a direct application of this model is precluded by a number of 
imperfections and irregularities present in real systems as illustrated in Fig. [T] These 
include a non-uniform waviness, mutual shift of individual yarns from layer to layer 
and most importantly a non-negligible porosity. Unlike finite element simulations 
which enable incorporating most of these imperfections directly through the for- 



mulation of a certain statistically equivalent periodic unit cell (Zeman and Sejnoha 



2004 



Sejnoha and Zeman, 2008), the MT method has only limited means which are, 
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Fig. 5. Ideal periodic unit cell: (a) cross-section, (b) three-dimensional view, (c) approx- 
imate distribution of inclination angles, (d) example of a real distribution of inclination 
angles 

nevertheless, still sufficient when quick estimates of the effective "bulk" response is 
needed. These will discussed next in the framework of the two-step homogenization 
procedure examined already in the previous section. 



4-1 Effective conductivities of plain weave textile composite ply 



Consider a simple plain weave fabric ply in the absence of porous ply. At this level, 
the carbon fiber tow is treated as a homogeneous phase with known material proper- 
ties bonded again to an isotropic carbon matrix. In order to address the infiuence of 



various geometrical fiaws, the approach proposed in (Skocek et al, 2008) is adopted. 
This involves: 



(1) determination of an ideal geometrical model to asses the volume fraction of the 
fiber tow 

(2) determination of the optimal shape of an equivalent ellipsoidal inclusion sub- 
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stituting the fiber tow in the solution of the Eshelby problem 
(3) proper evaluation of orientation dependent quantities from Eq. (16) to account 
for a non-uniform waviness along the fiber tow path 



4.1.1 Ideal geometrical model - Fig. ^b) 



The three-dimensional geometrical model adopted in the present work (Kuhn and 



Charalambides , 1999) is characterized by four parameters: the tow wavelength 2a, 
the tow height b, tow spacing g and the layer thickness h. To formulate one particular 
"ideal" representative, a tedious image analysis of a number of sections of a real 
textile ply such as the one in Fig. [l](ai) was carried out. The averages of the basic 
geometrical parameters presented in Table |3] were used to construct the required 
geometrical model. 
Table 3 



Quantification of PUCl parameters (Tomkova 2004) 



Statistics [/im] 


a 


h 


b 


9 


Average 


2250 


300 


150 


400 


Standard deviation 


155 


50 


20 


105 



4.1.2 Optimal shape of the equivalent ellipsoidal inclusion 



An extensive numerical study was performed in (Skocek et al, 2008) to conceive 
how the Mori-Tanaka predictions are influenced by a "random" deviation of basic 
geometrical parameters of real systems from their ideal representative introduced in 
the previous section. The results revealed a certain correlation between the model 
parameters and "optimal" shape of an equivalent ellipsoidal inclusion characterized 
by three semi-axes i^i, .^2) Ca- When setting ^1 = 1 (recall that the Eshelby solution 
depends only on the mutual ratio of the ellipsoid semi-axes), it was concluded that 
the ^2 parameter is strongly correlated with g/a ratio, while it is almost independent 
of b/a value. An analogous trend could be observed between ^3 and b/a parameter. 
This led to following semi-empirical formulas employed also in this study 

3(7 1 6 , , 

a 10 oa 

where the necessary parameters a, b, g are available from Table |3| 

Remind, however, that these relations were originally derived to match effective 
elastic properties. To either confirm or displace their validity in the solution of heat 
conduction problem we further assumed an equivalent inclusion in the form of an 
infinite elliptic cylinder with elliptical cross-section estimated directly from actual 
cross-section of tows in the representative model giving 

6 = 00, 6 = 12.3, ^3 = 1. (21) 
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4.1.3 Orientation averaging 



There are three distinct routes available in this study to reflect the variation of 
the inclination angle along the fiber tow path. If we consider directly the simplified 
geometrical model in Fig. IsFb), the joint probability density function g((f),9,P) re- 



sults from the harmonic shape of the centerline as described by Eq. (19). Applying 
the change of variable formula (Rektorys, 1994, Section 33.9), we obtain after some 
algebra the expression of the probability density in the form 

2a l + tan2(^) 



if 



otherwise 



0,/3 = and - a < 9 < a, 



where 



a 



arctan 



2a 



Assuming simply a uniform distribution of inclination angles the joint probability 
density function attains the form 



if 



2a 

otherwise. 



0, /5 = and - a < 9 < a, 



Both functions are plotted in Fig. |5](c) for comparison. In this study, the latter 
function was adopted for simplicity. Next, let D represents an orientation dependent 
quantity in Eq. (16) 



written, for the warp system, as 

((D™P)) 

and similarly for the fill system we get 



;D)) = ((X2A2))-Xi((A2)), 

g{<P,9,P)D{0,9,0)d9, 



9[ 



.9,(3)D{-,9,0)d9. 



(22) 
(23) 

(24) 



Finally, following ( Skocek et al. , 2008 ) , the resulting homogenized stiffness matrix 
given by Eq. (16) then becomes 



warp\ 



+ D 



(25) 



One may also suggest to model the plain weave fabric as a three-phase composite 
with warp and fill systems of yarns being considered as two distinct phases. The 
homogenized conductivity matrix then attains a slightly different form 



X 



warp A warp 



fill A fill 



Xi 



1 warp 



+ A 



fill 



(26) 
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where 



1 warp -pwarp 



cil + ^ (t^'^^'p + T^" 



1 -1 



-fill 



cil + ^ (t™p + T'^" 



1 -1 



-warp 



-fill 



T™^''P(0,^,0), 



-fill 



(27) 
(28) 



However, the differences in predictions provided by Eqs. (25) and (26) are, as seen 
in Table |4| neghgible. 

Improvements when compared to the assumed ideal geometry are contained in the 
third route which allows us to introduce the non-uniform waviness and to some 
extent also the mutual shift of individual layers by utilizing histograms of inclination 
angles shown in Fig. Istd). These are derived from centerlines of individual fiber tows 



described in detail in (Vopicka, 2004). Representing the joint probability density 
function by these histograms, the contribution to the effective conductivity matrix, 
e.g. from the warp direction Eq. (23), reads 



V warp \ 



5]p,D(O,0„O), 



(29) 



i=l 



where m denotes the number of sampling values. The discrete angles 6i and proba- 
bilities Pi follow directly from the image analysis data. Eleven such histograms as- 
sociated with several sections measured along individual plies were considered. The 
resulting averages together with the estimates provided by the simplified distribution 
functions are summarized in Table |4j For the solution of the Eshelby problem of an 
ellipsoidal inclusion in an isotropic matrix we refer the interested reader to (Hatta 



and Taya 


1986 


Jeong et al, 


1998 



4-2 Effective conductivities of homogenized porous matrix 



When carefully observing Fig. [T[a4) we identify three more or less periodically re- 
peating geometries further displayed in Fig. [6j These segments readily confirm the 
need for the proposed two step homogenization procedure as the ideal representative 
plotted in Figs. |6](a)(b) (already analyzed in the previous section) cannot be used to 
represent the entire composite. Instead, the second homogenization step is required 
to account for the presence of large vacuoles evident from Figs. [6]^c)-(f). 

Owing to the orthogonal arrangement of tows in the ideal (representative) model the 
new homogenized matrix employed in the second homogenization step is no longer 
isotropic. Thereby, the Eshelby solutions used so far are not directly applicable. 
Instead, the S tensor is found by imagining an equivalent ellipsoidal inclusion in 
an infinite matrix being orthotropic. The corresponding Laplace equation governing 
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(c) 



(f) 



Fig. 6. Homogenization on meso-scale: (a)-(b) PUCl representing carbon tow-carbon ma- 
trix composite, (c)-(d) PUC2 with vacuoles aligned with delamination cracks due to slip 
of textile plies, (e)-(f) PUC3 with extensive vacuoles representing the parts with textile 
reinforcement reduction due to bridging effect in the middle ply 



the steady state heat conduction problem is provided by 



Xi 



Introducing the following substitutions 



XlXl] X2 — y/X2X2] X- 



X3X3 



allows us to convert Eq. ( 30 ) into 



d^9 d^9 d^9 



which eventually leads to 

Sii = 



dXi 



dxl 



dxi 



(30) 



(31) 



(32) 



2^X1 X2 X3 Jo {{^iy/Xi + s) As 



ds, 



where 



As 



(33) 



(34) 



These equations then formally resemble those derived for the case of an isotropic 



in (Jeong et ai, 1998). 



matrix. The solution of Eq. (33) given in the form of elliptic integrals is available 



A simple example of an isotropic void (ellipsoidal inclusion with ^ = (1, 1.5, 2), = 
0.2) surrounded by an orthotropic matrix {xTi = 20, X22 = X33 = 2) was considered 
to acknowledge correctness of Eq. ( [33) ). Fig. [t] compares the MT predictions with 
the finite element results found for a hexagonal arrangement of voids under periodic 
boundary conditions (Tomkova et ai, 2008[ ). 
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Fig. 7. Case study: effective conductivities in [Wm -"^K ^] as a function of the volume 
fraction of the porous phase 

With these encouraging results at hand we proceeded with the analysis of real sys- 
tems. Unlike the micro-scale, the method of observation and measuring tools pro- 



vided by the LUCIA G software (LIM, 2008) were utilized here to approximate the 
shape of individual vacuoles in Figs. [6|^c)-(f). Since only two-dimensional (2D) im- 
ages were supplied, the voids were assumed to be well approximated by an oblate 
spheroid defined in Table |5] for both types of representative periodic unit cells in 
Figs. [6]^d)(f). To estimate the mesoscopic effective properties the voids were in- 
troduced into the homogenized matrix derived in the 1st homogenization step by 
combining Eqs. (25) and (29) (the last two columns in Table |4]). 



A word of caution, however, is appropriate when dealing with 2D images only. While 
the shape of the inclusion acquired from 2D images may play a minor role in final 
predictions of the effective properties, the volume fraction of a relevant heterogeneity 
also estimated from 2D images may prove much more important. This is documented 
in Table |4] listing the predicted effective properties of a plain weave fabric free of 
pores for two different volume fractions of the homogenized fiber tow. The volume 
fraction 0.53 corresponds to a representative 3D meso-structure in Fig. [5|b), whereas 
the value of 0.78 follows from the corresponding 2D image in Fig. [61 



Table 4 

Effective thermal conductivities of a plain weave fabric without porosity in [Wm~^K~^] 



Method 


Fiber tow 
volume 
fraction 


Thermal conductivity 


Eqs. (25)&(23) 


Eq. ^ 


Eqs. (25)&(29) 


warp/fill 


trans. 


warp/fill 


trans. 


warp/fill 


trans. 


MT 
Eq. (20) 


0.53 


8.18 


3.08 


7.93 


3.07 


8.14 


3.16 


0.78 


9.23 


2.37 


8.68 


2.34 


9.18 


2.48 


MT 
Eq. (21) 


0.53 


8.35 


2.94 


8.26 


2.92 


8.31 


3.02 


0.78 


9.35 


2.30 


9.15 


2.27 


9.29 


2.42 


3D FEM 


0.53 


8.13 


3.18 
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The mesoscopic effective conductivities derived for individual geometries in Fig. [6] 
are summarized in Table |5] for both volume fractions of the fiber tow. Note that 
the results corresponding to a representative model denoted as PUCl are essentially 
those stored in Table |4] which of course served as a point of departure for the 2nd 
homogenization step performed for models PUC2 and PUC3. 

Table 5 

Effective thermal conductivities of the porous textile plies and laminates in [Wm~^K~^] 



Fiber tow in 1st step 


Void 


Thermal conductivity 


vol. frac. 


shape 


PUC 


shape 


vol. frac. 


warp/fill 


transverse 


(histograms) 
0.53 




PUC2 


3;3;1 


0.07 


7.40 


2.75 


Eq. ( 


20 


) 


PUC3 


1.6; 1.6; 1 


0.15 


6.44 


2.51 




PUC2 


3;3;1 


0.07 


7.55 


2.63 


Eq. ( 


21 


) 


PUC3 


1.6; 1.6; 1 


0.15 


6.57 


2.40 


(histograms) 
0.78 




PUC2 


3;3;1 


0.07 


8.32 


2.18 


Eq. ( 


20 


) 


PUC3 


1.6; 1.6; 1 


0.15 


7.21 


1.99 




PUC2 


3;3;1 


0.07 


8.42 


2.13 


Eq. ( 


21 


) 


PUC3 


1.6; 1.6; 1 


0.15 


7.29 


1.95 



Table 6 

Effective thermal conductivities of the laminate [Wm'^K^-*^] 



Method 


Fiber tow in 1st step 


Thermal conductivity 


shape 


volume fraction 


warp/fill 


transverse 


MT 
(histograms) 


Eq. (20) 


0.53 


7.26 


2.76 


0.78 


8.15 


2.18 


MT 


Eq. (21) 


0.53 


7.40 


2.64 


(histograms) 




0.78 


8.25 


2.13 


2D FEM (Tomkova et al. 


2008 


) 




0.78 


8.67 


1.66 


Measured (^Kubicar et al!, 
■ 


2002 








10 


1.6 



5 Macro-scale 



The final, clearly the most simple, step requires construction of the homogeneous 
laminated plate. The stacking sequence of individual periodic unit cells complies 
with that observed for the actual composite sample (Tomkova et al, 2008, Fig. 2), 



see also Fig. |Tfa4) identifying the PUC1/PUC2/PUC3/PUC1 stacking sequence. 
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While in-plane conductivities (warp/fill directions) were found from a simple arit- 
metic rule of mixture, the out-of-plane (transverse) conductivity followed from the 
inverse (geometric) rule of mixture. The resulting effective thermal conductivities are 
available in Table [6] comparing the MT predictions and experimental measurements 
presented in (Kubicar et a/. 2002). The FEM results obtained from 2D simulations 
in (Tomkova et a/., 2008) are provided for additional comparison. Note that the 



highlighted (bolt font) values of thermal conductivities stored in Tables [2]-[6] follow 
from what we would call an optimal or the most appropriate approach. 



6 Discussion and conclusions 



In order to realistically model complex plain weave textile laminates with three- 
dimensional, generally non-uniform texture of the reinforcements and significant 
amount of porosity we advocate to consider at least three levels of hierarchy - the 
level of fiber tow, the level of yarns and the level of laminate. On each level different 
resolutions of microstructural details are considered on individual scales for the 
formulation of an adequate representative model. The desired macroscopic effective 
properties of the laminate are then estimated with regard to two basic objectives: 

• to reflect the three-dimensional character of the composite at all scales, 

• to predict the effective conductivities as efficiently as possible. 

Unlike computationally tedious and extensive 3D finite element simulations the 
Mori-Tanaka averaging scheme appears as a reasonable candidate to comply with 
both objectives. Not only the fully explicit format of this method but a simple exten- 
sion of the Eshelby problem, at least in the case of the solution of heat conduction 
problem, to generally orthotropic reference medium (homogenized composite free of 
pores in our particular case) favors this technique. 

In this study, the hierarchical character of the analysis is presented in a totally 
uncoupled format. Therefore, each level is treated entirely independently purely 
upscaling the results from a lower to higher scale for subsequent calculations. An 
"optimal" procedure, which attempts to accommodate various sources of imperfec- 
tions observed in real systems is accompanied by several modifications involving 
mainly the meso-scale. 



Based on our previous study of effective elastic properties (Zeman and Sejnoha 



2004 Skocek et al. 2008 ) it was expected that at this level the "best" estimates 



of the effective conductivities would follow from the application of Eq. (20) to de- 
termined the shape of an equivalent inclusion for the fiber tow representation and 
histograms of fiber inclination angle to proceed with the orientation averaging step. 
Comparing various modifications (different type of inclusion, ideal path of the fiber 
tow) suggests, perhaps even intuitively, almost negligible sensitivity of the solution 
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of the heat conduction problem to mutual interlacing of individual tows in com- 
parison with the solution of the elasticity problem. This is mainly attributed to a 
relative flatness of the reinforcing yarns in individual plies of the laminate. In view 
of this, one may even offer the possibility of estimating the effective mesoscopic 
conductivities by simply assuming two systems of perpendicular fiber tows with no 
interlacing, thereby completely avoiding the orientation averaging step. But bear in 
mind that such a "drastic" simplification can hardly be generalized and is certainly 
not acceptable in the case of elasticity. To discriminate between various approaches 
is therefore difficult. 



Comparison with experimental measurements is in principle twofold but also in- 
conclusive. On the one hand it clearly supports the use of the proposed uncoupled 
multi-scale approach and the two-step homogenization scheme on individual scales 
to arrive at the predictions of the effective macroscopic thermal conductivities. Fur- 
thermore, at least quantitatively, the Mori-Tanaka method proved its applicability 
in the solution of complicated textile composites. These remarks have already been 



put forward in (Skocek et al. 2008) with regard to the problem of effective elastic 



properties. To judge, on the other hand, the pertinence and reliability of the MT 
method solely by comparing the predicted and measured values, which may deny 
it, is certainly insufficient. While all deficiencies of the presented homogenization 
strategy were openly discussed, errors associated with experimental measurements 
were not mentioned and are not available. 



In summary, focusing only on the quantitative perspective, the Mori-Tanaka method 
combined with popular multi-scale homogenization approach is viable and presents 
a suitable and efficient alternative to periodic homogenization typically based on 
finite element analysis. 
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